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r--s ■ Angular dependent torque measurements were performed on single crystals of CeCoIng, a heavy 

fvj ' fermion superconductor [Tco = 2.3 K), in the temperature, T, range 1.9 K < T < 20 K and magnetic 

fields H up to 14 T. A large paramagnetic effect is found in the normal state. Torque measurements 

^\, in the mixed state were also performed. The torque curves show sharp hysteresis peaks at = 90° 

^f^ • {6 is the angle between H and the c— axis of the crystal), a result of intrinsic pinning of vortices. 

^ , The anisotropy 7 = \Jmclma in the mixed state was determined from the reversible part of the 

t^^ . vortex contribution to the torque signal using Kogan's model [Phys. Rev. B 38, 7049 (1988)]. The 

T-H ' anisotropy 7 decreases with increasing magnetic field and temperature. The fact that 7 is not a 

\ constant points towards a multiband scenario in this heavy fermion material. 

'O \ 
5_i . A. Introduction 

^ , 

^ , The heavy fermion compound CcCoIns forms in the HoCoGas tetragonal crystal structure with alternating layers 

3 ' of Celna and CoIn2. It is superconducting at 2.3 K, the highest superconducting transition temperature TcO yet 

^ , reported for a heavy fermion superconductor,^ Considerable progress has been made in determining the physical 

' , ' properties of this material. The superconductivity in this material is unconventional. The presence of a strong 

-H ' magnetic interaction between the 4f moments and itinerant electrons allows the possibility of nonphonon mediated 

Q I coupling between quasiparticlesr'"^ Angular dependent thermal conductivity measurements show dj.2_y2 symmetry 

C^ ' which implies that the anisotropic antiferromagnetic fluctuations play an important role in superconductivity^ Bel et 

'""' . al. reported the giant Nernst effect in the normal state of CeCoIns , which is comparable to high Tc superconductors 

I ' in the superconducting state. ^ Non-Fermi liquid behavior was observed in many aspectsi^ The unconventional 

h>. , superconductivity and the similarity to high Tc superconductors attract great interest to study this system. 

^\ • Measurements of de Haas- van Alphen oscillations in both the normal and mixed states have revealed the quasi two- 

0^ [ dimensional nature of the Fermi surface and the presence of a small number of electrons exhibiting 3D behaviorAS 

C^ Through point-contact spectroscopy measurement, Rourke et alpiS found that in CeCoIns there are two coexisting 

>! ^ order parameter components with amplitudes Ai = 0.95 ± 0.15 meV and A2 = 2.4 ± 0.3 meV, which indicate 

r^ . a highly unconventional pairing mechanism, possibly involving multiple bands. This is very similar to the case of 

^-^ ' MgB2, where two coexisting s-wave gaps were found by the same technique. ^^ Thermal conductivity and specific heat 

■"■^ — , measurements made by Tanatar et al. have revealed the presence of uncondensed electrons, which can be explained 

C^ ■ by an extreme multi-band scenario, with a d-wave superconducting gap on the heavy-electron sheets of the Fermi 

Ch ] surface and a negligible gap on the light, three-dimensional pocketsji^ 

I • The presence of multibands (gaps) together with the reported field dependence of the cyclotron masses^ point 

' O towards a possible temperature and/or field dependent anisotropy in the superconducting state, which according to 

b^ the standard anisotropic Ginzburg-Landau (GL) theory is given by 7 = ^/rn^Jm^ — H^c2 /^c2 — ^c/^a = ^a/^c (c 
fj , and a are crystallographic axes, and m,Hc2,^, and ^ are the effective mass, upper critical field, penetration depth, 
J^' ■ and coherence length, respectively). Specifically, in the multiband scenario proposed by Rourke et aliiS and Tanatar 
et al.fi2i different gaps may behave differently in magnetic field, which may lead to a field dependent 7. Also an 
anisotropic gap may result in a temperature dependent 7. Reports up to date give values of the anisotropy of 
^ ' CeCoIns in the range 1.5 to 2.47. For example, Petrovic et al. have reported an anisotropy of at least 2, as estimated 
from the ratio of the upper critical fields Hc2 along the c and a directionsii Measurements of Hc2(^) at 20 mK give 
a value for the anisotropy of about 2.47pi^ Magnetization measurements of the lower critical field for temperatures 
between 1.5 and 2.1 K give the ratio of the out-of-plane and in-plane penetration depth Xc/Xa ~ 2.3 and of the 
in-plane and out-of-plane coherence length £,a/£,c ~ 1-5, which gives an anisotropy of 2.3 and 1.5, respectively^ 

Magnetic torque is a sensitive tool for probing the anisotropy. It has been successfully applied to investigate the 
highly anisotropic high temperature superconductors and also the less anisotropic materials such as MgB2.^^'^^ But 
all these previous torque measurements were made on materials which have negligible paramagnetism. On the other 
hand, CeCoIns is a magnetic superconductor, so it may have large paramagnetism which cannot be ignored in the 
study of the mixed state. Here, we report torque measurements on single crystals of CeCoIus both in the normal state 
and the superconducting state. Our results show large paramagnetism in this material in the normal state. Therefore, 
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we assume that there are two contributions to the torque signal in the mixed state of CeCoIns single crystals: one 
coming from paramagnetism and the other one coming from vortices. We determined 7 from the reversible part of 
the vortex signal and found that 7 is not a constant, instead, it is field and temperature dependent. This provides 
evidence that the picture in this unconventional superconductor is not a simple single band scenario, supporting the 
conclusions of Rourke et al>i£ and Tanatar et alm^ 

B. Experimental Details 

High quality single crystals of CeCoIns were grown using a flux method. The superconducting transition temper- 
ature Tc, defined as the value of T at zero resistivity, is 2.3 K. The surfaces of the single crystals were etched in 
concentrated HCl for several hours and then thoroughly rinsed in ethanol in order to remove the indium present on 
the surface. The mass of the single crystal for which the data are shown is 0.75 mg. Angular dependent measure- 
ments of the magnetic torque experienced by the sample of magnetic moment M in an applied magnetic field H , were 
performed over a temperature range 1.9 K < T < 20 K and applied magnetic field range 0.1 T < iJ < 14 T using a 
piezoresistive torque magnetometer. In this technique, a piezoresistor measures the torsion, or twisting, of the torque 
lever about its symmetry axis as a result of the magnetic moment of the single crystal. The sample was rotated in 
the applied magnetic field between H \\ c-axis (6* — 0°) and H \\ a-axis {0 = 90°) and the torques Tine and Tj^ec were 
measured as a function of increasing and decreasing angle, respectively, under various temperature - field conditions. 

The contributions of the gravity and puck to the total torque signal were measured and subtracted from it. To 
measure the background torque due to gravity, we measured the torque signal at different temperatures in zero applied 
magnetic field with the sample mounted on the puck. The gravity torque is almost temperature independent and 
it is negligible at high applied magnetic fields. However, as the applied magnetic field decreases, the total torque 
signal becomes smaller and the effect of gravity becomes important, hence should be subtracted from the measured 
torque. To determine the contribution of the puck to the measured torque, we measured the torque without the single 
crystal on the puck at different magnetic fields and temperatures. The magnitude of the torque of the puck increases 
with increasing magnetic field. Also, the contribution of the puck to the measured torque is much larger than the 
contribution of the gravity. Therefore, the former contribution should always be subtracted from the total measured 
torque signal. 

C. Results and Discussion 

Previous torque studies of Tl2Ba2Cu06+i^^ and MgE^i^ti^ systems have shown that the normal state torque 
coming from paramagnetism is small compared with the fiux-vortex torque, therefore, one could neglect the former 
contribution to the total torque signal measured in the superconducting state. However, this is not the case for 
CeCoIns, for which the paramagnetic torque signal in the normal state is comparable with the total torque signal 
measured in the superconducting state, as shown below. Hence, one needs to subtract the former signal from the 
latter one in order to determine the torque due to vortices. This is similar to the case of the electron-doped high- 
Tc cuprate Ndi_2.Cea;Cu04, where a large paramagnetic contribution from Nd ions is discussed seperately from a 
superconducting contribution>i& Therefore, we first discuss the field and temperature dependence of the magnetic 
torque in the normal state and show how we subtract this contribution from the measured torque in the mixed state, 
and then we return to the discussion of the torque signal in the mixed state and to the determination of the field and 
temperature dependence of the bulk anisotropy 

All the torque curves measured in the normal state, some of which are shown in Fig 1, are perfectly sinusoidal and 
can be well fitted with 

Tp{T,H,0) = A{T,H) sm2e, (1) 

where A is a temperature and field dependent fitting parameter. Indeed, note the excellent fit of the data of Fig. 1 
with Eq. (1) (solid lines in the figure). The field dependence of A/H at 1.9, 6, 10 and 20 K is shown in Fig. 2. The 
solid lines are linear fits to the data, which show that A/H is linear in H with a negligible y-intercept and a slope 
which increases with decreasing T. So A is proportional to H^, i.e., 

A{T,H)^C{T)H^, (2) 

with C a temperature dependent fitting parameter. 



Next we show that the torque measured in the normal state and given by Eq. (1) is a result of the paramagnetism. 
Indeed, the torque of a sample of magnetic moment M placed in a magnetic field H is given by 

fp{T,H) = a1 X H. (3) 

The resultant magnetic moment M can always be decomposed into a component parallel Mn and one perpendicular 
M± to the a6— plane of the single crystal. With the magnetic field H making an angle 9 with the c— axis of the single 
crystal, Eq. (3) becomes: 

Tp{T,H,9) ^ [M\\H cose - Mj_H sme]k. (4) 

On the other hand, the experimental relationship of the torque, given by Eq. (1) with the fitting parameter A given 
by Eq. (2), becomes 

Tp{T,H,6) ^ A{T,H)sm2e ^ 2C{T)H'^ sinOcose. (5) 

Therefore, with C = (Ci - C2)/2, Eqs. (4) and (5) give 

Af|| ^CiHsme = XaH\i, 

where Ci = Xa, the a-axis susceptibility, and 

Mi_ ^C2H cose = XcH^, (6) 

where C2 = Xc, the c-axis susceptibility. This shows that the torque measured experimentally is of the form: 

rp(r, iJ, e) = ^^1—^H^ sin 26* (7) 

The fact that A/H = (xo - Xc)H/2 = {M\\ - Mj_)/2 [see Eqs. (5) and (7)] shows that A/H plotted in Fig. 2 reflects 
the anisotropy of the magnetic moments along the two crystallographic directions, a and c, while its linear field 
dependence shows that the magnetic moments are linear in H, hence the susceptibilities along these two directions 
are field independent. The temperature dependence of the magnetic moments is given by the temperature dependence 
of the parameter C. Therefore, Eq. (7) shows that the T, H, and 9 dependences of the torque measured in the normal 
state of CcCoIus reflect its paramagnetism and the anisotropy of its susceptibility Ax = Xa ^ Xc along the a and c 
directions. 

To further check the consistency of the data and to precisely determine the paramagnetic value of the torque, we 
also measured the magnetic moment M of the same single crystal of CeCoIus using a superconducting quantum 
interference device (SQUID) magnetometer. The magnetic moments measured at 4, 6, 10, 15, and 20 K are plotted 
in the main panel of Fig. 3 as a function of the applied magnetic field for both H \\ c-axis {e = 0°) and H \\ a-axis 
{e = 90°). The magnetic moments for both field orientations are linear in H with M|| < M± for all temperatures 
measured, consistent with the torque data of Fig. 2 and with Eq. (6). Note that the units for the magnetic moments 
of Figs. 2 and 3 are different. We change the units and compare the results of the two types of measurements. For 
example, the torque measured at 6 K and 5 T gives A/H = -1.42 x 10""^ Am~^ Since A/H = (M|| - M±)/2, 
AM = M|| — M± — —2.84 x 10~^ Am~^. The SQUID measurements give at the same temperature and applied 
magnetic field M± = 7.2 x lO"'' emu = 7.2 x 10""^ Am^^ and Af|| = 4.2 x 10^'^ emu = 4.2 x 10"'^ Am^^; hence, an 
anisotropy AM = —3.03 x 10~^ Am^^. Therefore, the values of the magnetic moments obtained in the two types of 
measurements are within 5% of each other, an error well within our experimental error. 

The Inset to Fig. 3 is a plot of the susceptibilities along the two directions, calculated from the slopes of M{H) of 
Fig. 3. Clearly x{T) shows anisotropy with respect to the field orientation. Note that the susceptibilities for both 
directions increase with decreasing temperature. The continuous increase of x{T) in the investigated temperature 
range may be related with the non- Fermi liquid behavior due to the proximity to the quantum critical field. ^^ These 
values of Xc and Xo are consistent with the ones reported by other groups.^ 

The above study of the magnetic torque in the normal state has shown that the contribution of paramagnetism 
to the torque signal is very large, it has a quadratic H dependence, and also a T dependence [see Eqs. (1) and (2)]. 
Therefore, to extract the vortex torque in the mixed state, one needs to account for this paramagnetic contribution 
and subtract the two background contributions from the measured torque. The gravity and puck contributions to the 
measured torque were determined and subtracted as explained in the Experimental Details Section. The resultant 
torque includes the paramagnetic Tj, and the vortex t^, contributions and is plotted in the inset to Fig. 4. We need to 



mention here that the vortex and paramagnetic torque contributions have opposite signs since the magnetic moment 
representing the vortex torque is diamagnetic. 

The Inset to Fig. 4 is the angular dependent torque measured in the mixed state at 1.9 K and 0.3 T in increasing 
and decreasing angles. Again, this torque includes Tp and Ty. Note that t(6) displays hysteresis. This hysteretic 
behavior is similar to the behavior in high Tc superconductors and is a result of intrinsic pinning.^*' The reversible 
component of the torque is determined as the average of the torques measured in increasing and decreasing angle; i.e., 

^rev — vdec i Tine) / ■^ yj) 

A plot of Tre^(0), obtained from Tinc{(^) and Tj^ecW data is shown in the main panel of Fig. 4. The reversible 
component of the torque reflects equilibrium states, hence it allows the determination of thermodynamic parameters. 
In the three-dimensional anisotropic London model in the mixed state, the vortex torque r„ is given by Kogan's 
modcli^ We assume that the paramagnetic contribution in the mixed state is given by Eq. (1). Therefore, 

Trev{0) =Tp + Ty =asm20+— -^ TffT ^"^ \ TT m\ ( (^) 

where a is a fitting parameter, V is the volume of the sample, /ig is the vacuum permeability, Xab is the penetration 
depth in the afc— plane, 7 = yjm^ma^ (.(&) = (sin^ + 7^ cos'^ ^)^^^, 77 is a numerical parameter of the order of unity, 

and iJ^2 is the upper critical field parallel to the c-axis [-ffj,2 (1.9 K) — 2.35 T]. We define /3 = (/)o-ffT^/(167r^oA^(,). 

To obtain the field dependence of the anisotropy 7, we fit the torque data with Eq. (9), with a, /?, and 7 as three 
fitting parameters. The solid line in the main panel of Fig. 4 is the fitting result for T = 1.9 K and H = 0.3 T. We 
need to mention here that the value of the fitting paramter a is 20% smaller than what we would expect from the 
extrapolation of the normal state paramangetic torque data. Also, a has an H dependence with an exponent of 2.3 
instead of 2. So, either there is extra contribution from other physics which has a weak field dependence, in addition 
to the paramagnetism contribution, or maybe the paramagnetic contribution becomes smaller in the mixed state of 
CeCoIns. Further experiments are needed to clarify this issue. 

Figure 5 is a plot of the field dependence of /3. The Inset is an enlarged plot of the low field region. Note that f3 
displays linear behavior up to a certain field with no y-intercept, then it deviates from linearity at iJ « 0.5 T, and 
it increases fast in the high field region. Since, on one hand A should be field independent, and on the other hand 
Eq. (9) is valid only for applied magnetic fields much smaller than the upper critical field, i.e. H « Hc2{T) for a 
given temperature, we assume that 0.5 T, the field at which (3{H) deviates from linearity, is the cutoff field Hcut for 
the applicability of the above theory. The slope of I3{H) in the linear H regime gives A (T = 1.9 K) = 787 nm. This 
value is larger than previous reports, which give Xab = 600 nm from measurements using a tunnel diode oscillator— 
and Xab = 330 nm from magnetization measurements.^"* 

Next, we fix A to the three values given above and fit the Trevid) data with only two fitting parameters, a and 7. 
The resultant field dependence of 7 is shown in Fig. 6 for the different A values. The parameter 7 is first decreasing 
with increasing H, reaches a minimum at H — 0.5 T, and then increases with further increasing field. As mentioned 
above, the cutoff field is 0.5 T. The data for H > Hcut are not reliable due to the failure of Kogan's theory in this H 
region. So we conclude that the anisotropy 7 decreases with increasing field. We note that this field dependence of 7 
in CeCoIns is opposite with the one for MgB2,-^- in which 7 increases with increasing field. We found that the value 
of 7 is very sensitive to the value of A, i.e., the larger the value of A, the larger the value of 7, with no effect however 
on its H dependence. 

To study the temperature dependence of 7, we performed torque measurements in the mixed state at 1.9 K, 1.95, 
2.00 K in an applied magnetic field of 0.3 T, and determined 7 from Eq. (9). Figure 7 is a composite plot of the 
temperature dependence of the anisotropy. The squares give "f(T) determined from torque measurements as discussed 
above [with A (1.9 K) = 600 nm, A (1.95 K) = 670 nm, and A (2 K) = 740 nm taken from Ref (25)], the sohd 

circles give 7(T) calculated from the ratio of H^c2 I -^cS taken from previous reports^ while the triangles give 7(T) 
determined from ^J pc{T) / pa{T) measured in zero field [see Pc{T) and Pa{T) in Inset to Fig. 7]. Note that the overall 
trend is a decrease of the anisotropy with increasing T, with a stronger dependence around rc(0.3 T) = 2.23 K. The 
values of 7 obtained from torque measurements are well within the range previously reported. The fact that the 
anisotropy depends both on temperature and magnetic field could explain its relatively wide range of values reported 
in the literature. 

The field and temperature dependence of the anisotropy implies the breakdown of the standard anisotropic GL 
theory, which assumes a single band anisotropic system with a temperature and field independent effective-mass 
anisotropy. In fact, a temperature and field dependent anisotropy was previously observed in NbSe2r^ LuNi2B2C,'^ 
and MgB2rSi2& the latter one having a value of 7 similar with CeCoIns. Very recently, Fletcher et al. reported 
measurements of the temperature-dependent anisotropics (7^ and 7^) of both the London penetration depth A and 



the upper critical field of MgB2 .^'- Their main result is that the anisotropics of the penetration depth and Hc2 in MgB2 
have opposite temperature dependences, but close to T^ they tend to a common value [see Fig. 4 of Ref (30)]. This 
result confirms nicely the theoretical calculation of anisotropy based on a two-gap scenarioiSP So, the temperature 
dependence of 7 in CeCoIns and MgB2 could have a similar origin. The field dependence of 7 could also be related 
with a multiband structure. In this scenario, the response to an applied magnetic field may be different in the case of 
the d-wave superconducting gap on the heavy-electron sheets of the Fermi surface and the negligible gap on the light, 
three-dimensional pockets. 

D. Summary 

Torque measurements were performed on CeCoIns single crystals in both the superconducting and normal states. 
Two contributions to the torque signal in the mixed state were identified: one from paramagnetism and the other 
one from the vortices. The torque curves show sharp hysteresis peaks aX 9 = 90° {0 is the angle between H and 
the c— axis of the crystal) when the measurements are done in clockwise and counterclockwise directions. This 
hysteresis is a result of the intrinsic pinning of vortices, a behavior very similar to high transition temperature 
cuprate superconductors. The temperature and magnetic field dependence of the anisotropy 7 is obtained from the 
reversible part of the vortex torque. We find that 7 decreases with increasing magnetic field and temperature. This 
result indicates the breakdown of the Ginzburg-Landau theory, which is based on a single band model and provides 
evidence for a multiband picture for CeCoIns, which is highly possible due to its complex Fermi topology^*S 
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I. FIGURE CAPTIONS 

Figure 1. (Color online). Angular dependence of the paramagnetic torque Tp measured in the normal state of 
CeCoIns at different temperatures T and applied magnetic field H values. The solid lines are fits of the data with 
Eq. (1). Inset: Schetch of the single crystal with the orientation of the magnetic field H and torque t with respect 
to the crystallographic axes. 

Figure 2. (Color online). Field H dependence of A/H, where A is the fitting parameter in Eq. (1). The solid lines 
are linear fits of the data. 

Figure 3. Plot of the magnetic moment M vs applied magnetic field H, with H \\ c-axis (solid symbols) and H \\ a— 
axis (open symbols), measured at 4, 6, 10, 15, and 20 K. Inset: Susceptibility x vs temperature T, measured with 
H II c-axis and H \\ a-axis. 

Figure 4. (Color online). Angular 6 dependence of the reversible torque Trev, measured in the mixed state of 
CeCoIn5 at a temperature T of 1.9 K and an applied magnetic field H of 0.3 T. The solid line is a fit of the data with 
Eq. (9). Inset: 9 dependence of the hysteretic torque r, measured in increasing and decreasing angle at the same T 
and H. 

Figure 5. Magnetic field H dependence of the fitting parameter (3. The solid line is a guide to the eye. Inset: 
Enlarged plot of the low field region of the data in the main panel. 

Figure 6. (Color online). Field H dependence of the anisotropy 7 measured at 1.9 K. 

Figure 7. (Color online). Composite plot of the temperature T dependence of the anisotropy 7. The circles show 
the results obtained from the upper critical field data [Ref. 26], the triangles are obtained from the resistivity data 
shown in the Inset of this figure, and the squares are from the present torque data measured in an applied magnetic 
field of 0.3 T. Inset: T dependence of the in-plane pa and out-of-plane pc resistivities measured in zero field. 
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